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Abstract 

Spacetime geometries dual to arbitrary fluid flows in strongly coupled N = 4 
super Yang Mills theory have recently been constructed perturbatively in the long 
wavelength limit. We demonstrate that these geometries all have regular event hori- 
zons, and determine the location of the horizon order by order in a boundary deriva- 
tive expansion. Intriguingly, the derivative expansion allows us to determine the 
location of the event horizon in the bulk as a local function of the fluid dynamical 
variables. We define a natural map from the boundary to the horizon using ingoing 
null geodesies. The area-form on spatial sections of the horizon can then be pulled 
back to the boundary to define a local entropy current for the dual field theory in the 
hydrodynamic limit. The area theorem of general relativity guarantees the positivity 
of the divergence of the entropy current thus constructed. 
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1 Introduction 



Over the last few years, a special class of strongly coupled ci-dimensional conformal field 
theories have been "solved" via the AdS/CFT duality. Quite remarkably, the solution to 
these theories is given by the equations of d + 1 dimensional gravity (interacting with other 
fields) in AdS^+i spacetime. Since the long distance dynamics of any genuinely interacting 
field theory is well described by the equations of relativistic hydrodynamics, it follows as a 
prediction of the AdS/CFT correspondence that at long distances, the equations of gravity 
in an AdS^+i background should reduce to the (relativistic) Navier-Stokes equations in d 
dimensions. There is now substantial direct evidence for the connection between the long 
distance equations of gravity on AdS^+i spacetime and d dimensional relativistic fluid 
dynamics; c/., 0, |, |, § g |, 0, §, |, 0, @ §, § 0, §0§0 
g§g§g§gggggggg,gg, 0, 0, |3J, [3§ [37|, [3§ for a sampling of the literature on 



the subject. 

In particular, it was noted in |H| that the equations of pure gravity with a negative 
cosmological constant form a universal subsector in any theory of gravity on AdS spacetime. 
Following up on earlier work [^, [H], |23|, it was demonstrated in p3 (for AdSs) and 
more recently in |37j (for AdS4) that Einstein's equations in this universal sector may be 



recast, order by order in a boundary derivative expansion, into equations of motion for 
two collective fields, namely - the 'temperature' and the 'velocity'. These new equations 
of motion turn out to be simply the relativistic Navier-Stokes equations of fluid dynamics. 



The gravitational solutions of [|33] and [0 constitute an explicit map from the space 
of solutions of the hydrodynamic equations to the space of long wavelength gravitational 
solutions (which are asymptotically AdS).[] Subject to a regularity condition that we will 
discuss further below, the solutions of |33|, [37] are locally exhaustive in solution space i.e., 
all long wavelength solutions to Einstein's equations that lie nearby in solution space to a 
metric dual to a particular fluid flow are themselves metrics dual to slightly perturbed fluid 
flows. This at first sight surprising result is a consequence of the requirement of regularity. 



■""By 'long wavelength' solutions, we mean solutions whose spacetime variations are slow on a scale set 
by their respective boundary extrinsic curvature. Via the AdS/CFT dictionary this is the same as the 
requirement that the solutions vary slowly on the scale of the inverse temperature associated with the local 
energy density of the solution. 
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This requirement cuts down the 9-parameter space of Fefferman-Graham type solutions of 
AdS5 spacetime - parameterized by a traceless boundary stress tensor - to the 4-parameter 
set of solutions of fluid dynamics. 

We believe the local exhaustiveness of the gravity solutions dual to fluid dynamics, 
described in the previous paragraph, in fact generalizes to a global statement. We think it 
likely, in other words, that the solutions of [^, |37| in fact constitute all long wavelength 
asymptotically AdS solutions of gravity with a cosmological constant; we pause here to 
explain why. Every state in a conformal field theory has an associated local energy density 
and a consequent associated mean free path length scale l m f p , the inverse of the tempera- 
ture that is thermodynamically associated with this energy density. As a consequence of 
interactions every state is expected to evolve rapidly - on the time scale l m f P - towards 
local thermodynamical equilibrium, in an appropriate coarse grained sense,f] at the local 
value of the temperature. This approach to local equilibrium is not long wavelength in 
time and is not well described by fluid dynamics. The dual bulk description of this (short 
wavelength) phenomenon is presumably gravitational collapse into a black hole. On the 
other hand, once local equilibrium has been achieved {i.e., a black hole has been formed) 
the system (if un-forced) slowly relaxes towards global equilibrium. This relaxation process 
happens on length and time scales that are both large compared to the inverse local tem- 
perature, and is well described by fluid dynamics and therefore by the solutions of |33|. j37| . 
In other words it seems plausible that all field theory evolutions that are long wavelength in 
time as well as space are locally equilibriated, and so are well described by fluid dynamics. 
The discussions of this paragraph, coupled with the AdS / CFT correspondence, motivate 



the conjecture that the solutions of [|33| , p7| are the most general regular long wavelength 
solutions to Einstein's equations in a spacetime with negative cosmological constant in five 
and four spacetime dimensions respectively. 

We pause here to note two aspects of the solutions of [[33|. [37| that we will have occasion 
to use below. First, it is possible to foliate these solutions into a collection of tubes, each of 
which is centered about a radial ingoing null geodesic emanating from the AdS boundary. 
This is sketched in Fig. |] for a uniform black brane, where we indicate the tubes on a local 
portion of the spacetime Penrose diagram.^ As we will explain below, the congruence of 
null geodesies (around which each of our tubes is centered) yields a natural map from the 
boundary of AdS space to the horizon of our solutions. When the width of these tubes in 
the boundary directions is small relative to the scale of variation of the dual hydrodynamic 
configurations, the restriction of the solution to any one tube is well-approximated tube- 



2 The N — » oo limit of the field theory (dual to the classical limit in gravity) justifies this coarse graining 
and supresscs consequent fluctuations (which arc dual to quantum fluctuations in gravity). 

3 This is a causal diagram which captures the entire globally extended spacetime (note that in order 
for the boundaries to be drawn straight, the singularities are curved, as discussed in [p9[). For a realistic 
collapse scenario, described by the nonuniform solutions of (33j, only the right asymptotic region and the 
future horizon and singularity are present. 
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wise by the metric of a uniform brane with the local value of temperature and velocity. This 
feature of the solutions - the fact that they are tube-wise indistinguishable from uniform 
black brane solutions - is dual to the fact that the Navier-Stokes equations describe the 
dynamics of locally equilibrated lumps of fluid. 

Second, the gravitational solutions constructed in |33[] are regular everywhere away from 
a spacelike surface, and moreover the authors conjectured that this singularity is shielded 
from the boundary of AdS space by a regular event horizon. We will prove this conjecture 
by explicitly constructing the event horizon of the solutions of |33| order by order in the 
derivative expansion. It should be possible to carry out a parallel study for the solutions 



presented in |37[ for four dimensions. We will not carry out such a study here; however, 
aspects of our discussion are not specific to AdS 5 and can be used to infer the desired 
features of 2 + 1 dimensional hydrodynamics. We expect that the results of such a study 
would be similar to those presented in this paper. 

As we have explained above, we study the causal properties - in particular, the structure 
of the event horizon for the solutions presented in |33| 
various aspects of the dynamics - specifically, the entropy production 



We then proceed to investigate 
at this event 



horizon. In the rest of the introduction, we will describe the contents of this paper in some 
detail, summarizing the salient points. 





Fig. 1: Penrose diagram of the uniform black brane and the causal structure of the spacetimes 
dual to fluid mechanics illustrating the tube structure. The dashed line in the second figure 
denotes the future event horizon, while the shaded tube indicates the region of spacetime 
over which the solution is well approximated by a tube of the uniform black brane. 

As we have discussed above, provides a map from the space of solutions of fluid 
dynamics to a spacetime that solves Einstein's equations. The geometry we obtain out of 
this map depends on the specific solution of fluid dynamics we input. In this paper we 
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restrict attention to fluid dynamical configurations that approach uniform homogeneous 
flow at some fixed velocity and temperature at spatial infinity. It seems intuitively 
clear from the dissipative nature of the Navier-Stokes equations that the late time behaviour 
of all fluid flows with these boundary conditions will eventually become u^{x) = u$ 
and T(x) = T^; we assume this in what follows.^ The gravitational dual to such an 
equilibrated fluid flow is simply the metric of a uniformly boosted black brane. 

The causal structure of the uniform black brane is given by the Penrose diagram plotted 
in Fig.[l| (see In particular, the equation for the event horizon of a uniform black 



brane is well known. The event horizon of the metric dual to the full fluid flow is simply 
the unique null hypersurface that joins with this late time event horizon in the asymptotic 
future.Q It turns out to be surprisingly simple to construct this hypersurface order by 



order in the boundary derivative expansion used in ||33|| . In this paper, we perform this 
construction up to second order in derivatives. Within the derivative expansion it turns out 
that the radial location of the event horizon is determined locally by values and derivatives 
of fluid dynamical velocity and temperature at the corresponding boundary point. This is 
achieved using the boundary to horizon map generated by the congruence of ingoing null 
geodesies described above (see Fig. [[]).[] 

However, while locality is manifest in a derivative expansion, upon summing all orders 
we expect this local behaviour to transmute into a limited nonlocality: the radial position 
of the event horizon at a given point should be determined by the values of fluid dynamical 
variables in a patch of size 1/T centered around the associated boundary point. The 
ateleological behaviour of the event horizon is a surprising feature of these solutions and 
implies that the event horizon behaves as a 'membrane' whose vibrations are a local mirror 
of fluid dynamics. Our explicit construction of the event horizon of the metrics dual to 
fluid dynamics is one of the key results of our paper; cf., ( ^4|) . 

We now turn to a description of our second main result; the construction of an entropy 
current with non-negative divergence for a class of asymptotically AdS solutions of grav- 



4 This is true for conformal fluid dynamics in d spacetime dimensions for d > 3. Conformal fluid 
dynamics in 1 + 1 dimensions is not dissipative (for instance it is non-viscous since here the shear tensor 
does not exist). More generally, there are as many degrees of freedom in a traceless 1 + 1 dimensional stress 
tensor as in temperature and velocity fields in 1 + 1 dimensions. Consequently, the most general solution 
to 1 + 1 dimensional 'fluid dynamics' is simply given by T ++ = /(cr + ) and T__ = h(cr~) for arbitrary 
functions g and h. This solution describes left and right moving waves that maintain their shape forever, 
propagating non-dissipativcly at the speed of light. We thank A. Strominger discussions on this point. 

5 Note that for a generic hydrodynamic solution, the bulk spacetime has no manifest isometries; the 
event horizon is therefore not a Killing horizon. 

6 This map may be motivated as follows. Consider perturbing the fluid at a boundary point x^, e.g., 
by turning on some local operator of Yang Mills theory. This perturbation instanteneously alters all fluid 
quantities, including the entropy, at x^ . However, it only alters the geometry near the horizon at and 
within the lightcone emanating from x^ at the boundary. It is therefore plausible that local properties of 
the spacetime in the neighbourhood of ingoing geodesies that emanate from x^ capture properties of the 
fluid at x^. 
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ity, and its explicit evaluation for the solutions of [33] at second order in the derivative 
expansion.[] As we will see in § || it is possible to define a natural area (d — l)-form on any 
event horizon in a d+ 1 dimensional spacetime in general relativity. This form is defined to 
ensure that its integral over any co-dimension one spatial slice of the horizon is simply the 
area of that submanifold. It follows almost immediately from the definition of this form 
and the classic area increase theorems of general relativity that the exterior derivative (on 
the event horizon) of this (d— l)-form, viewed of as a top dimensional form on the horizon, 
is 'positive' (we explain what we mean by the positivity of a top form on the horizon in 

The positivity of the exterior derivative of the area (d — l)-form is a formally elegant 
restatement of the area increase theorem of general relativity that is local on the horizon. 
Hence we would like to link this statement to the positivity of the entropy production in 
the boundary theory. However, at least naively, the CFT fluid dual to our solutions lives 
at the boundary of AdS space rather than on its horizon. If we wish to study the interplay 
between the local notion of entropy of the fluid and the fluid equations of motion, it is 
important for these quantities to be defined on the same space. In order to achieve this, 
in § [| we use a congruence of null geodesies described above to provide a 'natural' map 
from the boundary to the horizon for a class of asymptotically AdS solutions of gravity 
(which include but are slightly more general than those of |33]). The pullback of the area 
(d — l)-form under this map now lives at the boundary, and also has a 'positive' exterior 
derivative. Consequently, the 'entropy current', defined as the boundary Hodge dual to the 
pull-back of the area (d — l)-form on the boundary (with appropriate factors of Newton's 
constant), has non- negative divergence, and so satisfies a crucial physical requirement for 
an entropy current of fluid dynamics. 

In § we then proceed to implement the construction described in the previous para- 
graph for the solutions of p3| . This enables us to derive an expression for the entropy 
current, Jg, with non-negative divergence, valid up to second order in the derivative ex- 
pansion. As a check of our final result, we use the equations of fluid dynamics to inde- 
pendently verify the non-negativity of divergence of our entropy current at third order in 
the derivative expansion. An example of an entropy current for a conformal fluid with 
non- negative divergence was first described in ||40|| . 

We also take the opportunity to extend and complete the analysis presented in [JEJ to 
find the most general Weyl covariant two derivative entropy current consistent with the 
second law. Note that the requirement of pointwise non-negativity of the entropy produc- 
tion - which we impose as a physical constraint of acceptable entropy currents - carries 
useful information even within the derivative expansion, though this is a little subtle to 
unravel. In particular, in § M we present a parameterization of the most general (7 param- 



7 We have been informed that A. Strominger and S. Hartnoll have independently discussed the con- 
struction of a positive divergence entropy current on the horizon, utilizing Brown York stress tensor. 
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eter) class of Weyl invariant candidate entropy currents that has the correct equilibrium 
limit, to second order in the derivative expansion. We also demonstrate that only a five 
dimensional sub-class of these currents is consistent with the requirement of pointwise 
non- negativity of d^Jg to third order in derivatives. We then turn our attention to the 
arbitrariness of our gravitational construction of the entropy current and demonstrate that 
there appears to be a two parameter family of physically acceptable generalizations of this 
bulk construction (associated with physically acceptable generalizations of the boundary 
to horizon map and also the generalisations of the area (d — l)-form itself). As a result, 
we conclude that the gravitational construction presented in this paper yields a two di- 
mensional sub-class in the five dimensional space of entropy currents with non-negative 
divergence. It would interesting to understand directly from field theory, what principle 
(if any) underlies the selection of this distinguished class of entropy currents. It would 
also be interesting to investigate whether the remaining positive entropy currents may be 
obtained from a generalized gravitational procedure, perhaps involving apparent and other 
local horizons.0 

This paper is organized as follows. In § [| below, we present an order by order construc- 
tion of the event horizon in a class of metrics that include those of |33[]. Following that, in 
§ 0, we present our construction of a local entropy (d — l)-form for the event horizons in 
gravity, and then implement our constructions in detail for the class of metrics studied in 
§ 0. In § |5|, we specialize the results of § [2] and § to the metrics dual to fluid dynamics 
p3| , using a map for translating horizon information to the boundary developed in § |j. We 
obtain explicit formulae, to second order in the derivative expansion, for the event horizon 
and the entropy current in the geometries of [33 . In § ^, we explain in detail the nature 
of the constraint imposed on second order terms in the expansion of the entropy current 
by the requirement of non-negativity of entropy production at third order in the derivative 
expansion. We also investigate the relationship between the geometrically constructed en- 
tropy current and the general entropy current of non-negative divergence generalizing the 
analysis of [f40|. Finally, in § ^, we end with a discussion of our results and open questions. 
Some technical results regarding the computations are collected in various Appendices. 



2 The Local Event Horizon 

As we have explained in the introduction, in this paper we will study the event horizon of 
the metrics dual to fluid dynamics presented in |H| . In that reference the authors construct 
an explicit classical spacetime dual to an arbitrary solution of fluid dynamics, accurate to 
second order in the boundary derivative expansion. While the explicit solutions of [|33| are 



rather involved, we will see below that the structure of the event horizons of these solutions 



3 We thank A. Strominger and M. Van Raamsdonk for discussions on this point. 
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are insensitive to many of these details. Consequently, in this section we will describe the 
metric of [33| only in general structural form, and carry out all our computations for an 
arbitrary spacetime of this form. In § |5] we will specialize these calculations to the detailed 
metrics of p3|. We start by presenting a geometric interpretation for the coordinate system 



used in [R3 



2.1 Coordinates adapted to a null geodesic congruence 

Consider a null geodesic congruence (i.e., a family of null geodesies with exactly one 
geodesic passing through each point) in some region of an arbitrary spacetime. Let E 
be a hypersurface that intersects each geodesic once. Let x^ be coordinates on E. Now 
ascribe coordinates (p, x M ) to the point at an affine parameter distance p from E, along 
the geodesic through the point on E with coordinates x M . Hence the geodesies in the 
congruence are lines of constant x M . In this chart, this metric takes the form 

ds 2 = -2u fl (x)dpdx fi + Xtiu{p, x)dx fl dx u , (2.1) 

where the geodesic equation implies that is independent of p. It is convenient to gener- 
alize slightly to allow for non-affine parametrization of the geodesies: let r be a parameter 
related to p by dp/dr = S(r,x). Then, in coordinates (r, x), the metric takes the formf] 

ds 2 = Gmn dX M dX N = -2 U(l {x) S{r, x) dr dx* + x^O, x) dx^ dx v (2.2) 

Note that E could be spacelike, timelike, or null. We shall take E to be timelike. 

This metric has determinant —S 2 x fll/ u fl u u det Xi where x^ v is the inverse of Xp,v- Hence 
the metric and its inverse will be smooth if S, and x^v are smooth, with S ^ 0, x^u 
invertible, and x^ v u n timelike. These conditions are satisfied on, and outside, the horizons 
of the solutions that we shall discuss below. 



2.2 Spacetime dual to hydrodynamics 

The bulk metric of |33|] was obtained in a coordinate system of the form ( |2.2j ) just described, 
where the role of E is played by the conformal boundary and the null geodesies are future- 
directed and ingoing at the boundary. The key assumption used to derive the solution 
is that the metric is a slowly varying function of x M or, more precisely, that the metric 

9 We use upper case Latin indices {M, N, ■ ■ ■ } to denote bulk directions, while lower case Greek indices 
{(i, ■ ■ ■ } will refer to field theory or boundary directions. Furthermore, we use lower case Latin indices 
{a, b, ■ • • } to denote the spatial directions in the boundary. Finally, we use (x) to indicate the depen- 
dence on the four coordinates . Details regarding the conventions used in this paper can be found in 
Appendix 
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functions have a perturbative expansion (with a small parameter e): 



S(r,*) = l-£e fe # , (2.3) 



oo 

X ,v(r, x) = -r 2 f(b r) u u +r 2 P„+^ e k ( S <*> r 2 P„ v + sf ] u v + u, + jj? u v + tjg) 

k=l 

(2.4) 

The function /(?/) above has the form / = 1 — ^; however, the only property of / that we 

will use is that /(l) = 0. The remaining functions (si fe \ s£ . . .) are all local functions of 
the inverse temperature b(x) and the velocity u^{x) and the coordinate r, whose form was 
determined in [33| and is indicated below in ( |5.1| ) and given explicitly in Appendix |A|; we 



however will not need the specific form of these functions for the present discussion. As 
far as the calculations in this section are concerned, the expressions s£\ sf*\ Sc , jji and 
t$ may be thought of as arbitrary functions of r and x^. The tensor P^ v = rj^ + u v is 
a co-moving spatial projector. 

In the above formulae, e is a formal derivative counting parameter. Any expression 
that multiplies e k in ( j2.3|) and Q2.4j) is of k th order in boundary field theory derivatives. 
Note that any boundary derivative of any of the functions above is always accompanied 
by an additional explicit power of e. As in |]33f| , all calculations in this paper will be 
performed order by order in e which is then set to unity in the final results. This is a good 
approximation when field theory derivatives are small in units of the local temperature. 

As we have explained in the Introduction, the metrics presented in [Q simplify to 
the uniform black brane metric at late times. This metric describes a fluid configuration 
with constant and b. As the derivative counting parameter e vanishes on constant 
configurations, all terms in the summation in (|2.3|) and ( |2.4j) vanish on the uniform black 



brane configuration. The event horizon of this simplified metric is very easy to determine; 
it is simply the surface r = \. Consequently, the event horizon 7i of the metric ( |2.2j ) has 
a simple mathematical characterization; it is the unique null hypersurface that reduces 
exactly, at infinite time to r = -r. 

In § ^]3|we will describe a local construction of a null hypersurface in the metric ( |2.2| ). 
Our hypersurface will have the property that it reduces exactly to r = 1/6 when w M and b 
are constants, and therefore may be identified with the event horizon for spacetimes of the 
form (|2.2|) that settle down to constant and b at late times, as we expect for metrics 
dual to fluid dynamics. We will evaluate our result for the metrics of |H§ in § || where we 
will use the explicit expressions for the functions appearing in fl2.2|) . 
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2.3 The event horizon in the derivative expansion 

When e is set to zero and b and w M are constants, the surface r = r is a null hypersurface in 
metrics ( |2.2|) . We will now determine the corrected equation for this null hypersurface at 
small e, order by order in the e expansion. As we have explained above, this hypersurface 
will be physically interpreted as the event horizon TC of the metrics presented in . 

The procedure can be illustrated with a simpler example. Consider the Vaidya space- 
time, describing a spherically symmetric black hole with ingoing null matter: 

ds 2 = -(l-'^^-) dv 2 + 2dvdr + r 2 dQ 2 . (2.5) 



Spherical symmetry implies that the horizon is at r = r(v), with normal n = dr — rdv. 
Demanding that this be null gives r(v) = 2m(v) + 2r(v) r(v), a first order ODE for r(v). 
Solving this determines the position of the horizon non-locally in terms of m(v). However, 
if we assume that m(v) is slowly varying and approaches a constant for large v, i.e., 

m(v) = O(e) ,mm = 0(e 2 ), etc., and lim m{v) = uiq (2-6) 

v — >oo 

then we can solve by expanding in derivatives. Consider the ansatz, r = 2m + amm + 
b m m 2 + c m 2 m + . . . , for some constants a, b, c, . . . ; it is easy to show that the solution for 
the horizon is given by a = 8, b = 64, c = 32, etc.. Hence we can obtain a local expression 
for the location of the horizon in a derivative expansion. 

Returning to the spacetime of [Q, let us suppose that the null hypersurface that we 



are after is given by the equation 

Sn(r,x) = , with Sn(r,x) = r — r#(x) . (2.7) 
As we are working in a derivative expansion we take 

M*) = ^y + 5> fc r (fc) (x) (2-8) 

Let us denote the normal vector to the event horizon by by definition, 

i A = G AB d B S H {r,x) (2.9) 

which also has an e expansion. We will now determine V(k)(x) and ^)(^ M ) order by order 
in e. In order to compute the unknown functions rn^(x) we require the normal vector £ A 
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to be null, which amounts to simply solving the equation 

G AB d A S n d B S n = (2.10) 
order by order in perturbation theory. Note that 

oo 

dS n = dr - edf.rH dx^ where e d^r H = - — d^b + ^ e " +1 <V (n) • (2.11) 

n=l 



In particular, to order e n , only the functions 77 m ) for m < n — 1 appear in ( 2.11 ). However, 
the LHS of ( |2.10| ) includes a contribution of two factors of dr contracted with the metric. 
This contribution is equal to G rr evaluated at the horizon. Expanding this term to order 
e n we find a contribution 



1 

r (n) 



K\ b 

where n\ is defined in ( [2.1 5|) below, together with several terms that depend on 77 m ) for 



Til 



< n — 1 . It follows that the expansion of Q2.10D to n th order in e yields a simple algebraic 



expression for r( n ), in terms of the functions rm, r( 2 ), • ■ ■ , ?"(n-i) which are determined from 
lower order computations. 

More explicitly, equation (|2.10|) gives G rr — 2ed il r H G r ^ / + e 2 d^r H d u r H G^ = 0, with 
the inverse metric G MN given by: 

-S 2 u tl u u x^ u ' S 2 u il u v x liv ' -S 2 u^x^ 

(2.12) 

where the 'inverse d-metric' is defined via = Hence the expression for the 

location of the event horizon (|2.10|) to arbitrary order in e is obtained by expanding 

= ^ (l - 2 e S X a % d a r H - e 2 S 2 ( X Q/3 X lS ~ X° 7 X P& ) « 7 u 5 d a r H d p r H ) 

(2.13) 

to the requisite order in e, using the expansion of the individual quantities S and r# 
specified above, as well as of x^ v ■ 
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2.4 The event horizon at second order in derivatives 



The equation Q2.1C ) is automatically obeyed at order e°. At first order in e we find that 
the location of the event horizon is given by r = withF] 

rg ) (i) = ^ + r ( i ) (x) = i + « 1 (sP-lrfd^ . (2.14) 



where we define 



At next order, 0(e 2 ), we find 



1 d n 



K m dr 1 



(r 2 f(br)) 



(2.15) 



r g)( x ) = I + Kl I s « + drsPrg - | (1 - s^) u" d,b + + 2 9 M r (1) 



2 



where we haveQ 

= u^u" + rT and if = diag(-l, 1, 1, 1) 



(2.16) 



As all functions and derivatives in (|2.14j) and ( [2. 16 ) are evaluated at r = 1/6 and the point 



x M and we retain terms to O (e) and O (e 2 ) respectively. 

It is now simple in principle to plug ( |2.16| ) into ( |2.2| ) to obtain an explicit expression 
for the metric H^ v of the event horizon.^ We will choose to use the coordinates x M to 
parameterize the event horizon. The normal vector £ is a vector in the tangent space of 
the event horizon (this follows since the hypersurface is null), i.e., 

£ A tAt = n^^- + n r ^- , (2.17) 
^ dX A dxf dr ' V ; 



10 We have used here the fact that = and tffi = which follow from the solution of J3j| . We 

also restrict to solutions which are asymptotically AdSs in this section. 

n It is important to note that in our expressions involving the boundary derivatives we raise and lower 
indices using the boundary metric ?7 M „; in particular, = rj^ u u v and with this defmtion = — 1. 



12 There are thus three metrics in play; the bulk metric defined in (2.2), the boundary metric which is 
fixed and chosen to be rj^ v and finally the metric on the horizon Ti, H^ Vl which we do not explicitly write 
down. As a result there are differing and often conflicting notions of covariance; we have chosen to write 
various quantities consistently with boundary covariance since at the end of the day we are interested in 
the boundary entropy current. 
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which is easily obtained by using the definition ( |2.9| ) and the induced metric on the event 
horizon; namely 



" = (i + «?> + (4 1} ) 2 + 4 2) ) «" - ^(* (1) r (i« + ^ 

+ ^ (i + «?> - »?>) + ^ (i - + i< 2) - ^r (1) 



:2.i8i 



Before proceeding to analyze the entropy current associated with the local area-form 
on this event horizon, let us pause to consider the expression ( 2.16|) . First of all, we see 



that for generic fluids with varying temperature and velocity, the radial coordinate r = ru 
of the horizon varies with x 11 , which, to the first order in the derivative expansion, is given 
simply by the local temperature. The constraints on this variation are inherited from the 
equations of relativistic fluid dynamics which govern the behaviour of these temperature 
and velocity fields, as discussed above. Note that the variation of r# at a given x % and as a 
function of time, can of course be non-monotonic. As we will see in the next section, only 
the local area needs to increase. This is dual to the fact that while a local fluid element 
may warm up or cool down in response to interacting with the neighbouring fluid, the local 
entropy production is always positive. An example of the behaviour of rn{x) is sketched 
in the spacetime diagram of Fig.0, with time plotted vertically and the radial coordinate 
as well as one of the spatial x % coordinates plotted horizontally. 



3 The Local Entropy Current 

Having determined the location of the event horizon, it is a simple matter to compute 
the area of the event horizon to obtain the area of the black brane. However, as we wish 
to talk about the spatio-temporal variation of the entropy, we will first describe entropy 
production in a local setting. This will allow us to derive an expression for the boundary 
entropy current in § 

3.1 Abstract construction of the area (d — l)-form 

In this brief subsection we present the construction of the area d — 1 form on the spatial 
section of any event horizon of a d + 1 dimensional solution of general relativity. 

First, recall that the event horizon is a co-dimension one null submanifold of the d+ 1 
dimensional spacetime. As a result its normal vector lies in its tangent space. The horizon 
generators coincide with the integral curves of this normal vector field, which are in fact null 
geodesies^ that are entirely contained within the event horizon. Let us choose coordinates 

13 This follows from the fact that the event horizon is the boundary of the past of future infinity T + 
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Fig. 2: The event horizon r = rn{x^) sketched as a function of the time t and one of the spatial 
coordinates x (the other two spatial coordinates are suppressed). 

(X,a a ), with a = 1, •■■ ,d— 1, on the event horizon such that a a are constant along 
these null geodesies and A is a future directed parameter (not necessarily affine) along the 
geodesies. As d\ is orthogonal to every other tangent vector on the manifold including 
itself, it follows that the metric restricted on the event horizon takes the form 

together with the fact that boundaries of causal sets are generated by null geodesies Q. We pause here 
to note a technical point regarding the behaviour of the horizon generators: While by definition these null 
geodesies generating the event horizon have no future endpoints [ ft2f , they do not necessarily remain on the 
event horizon when extended into the past. This is because in general dynamical context, these geodesies 
will have non-zero expansion, and by Raychaudhuri's equation they must therefore caustic in finite affine 
parameter when extended into the past. Hence, although the spacetime, and therefore the event horizon, 
are smooth, the horizon generators enter the horizon at points of caustic. However, since the caustic locus 
forms a set of measure zero on the horizon, in the following discussion we will neglect this subtlety. 
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ds 2 = g a bda a da b (3.1) 

Let g represent the determinant of the (d — 1) x (d — 1) metric g a b- We define the entropy 
(d — l)-form as the appropriately normalized area form on the spatial sections of the 
horizon^ 

a = — — y/gda 1 A da 2 A ... A da d ~ x (3.2) 

The area increase theorems of general relativity^ are tantamount to the monotonicity of 
the function g, i.e., 

I £ < 3 ' 3 > 

which of course leads to 

da = dX Ada 1 Ada 2 ... A da^ 1 > . (3.4) 

We have chosen here an orientation on the horizon H, by declaring a d-form to be positive 
if it is a positive multiple of the d-form dX A da 1 A da 2 ... A <ia d_1 . 

3.2 Entropy (d — l)-form in global coordinates 

The entropy (d — l)-form described above was presented in a special set of a a coordinates 
which are well adapted to the horizon. We will now evaluate this expression in terms of a 
more general set of coordinates. Consider a set of coordinates x^ for the spacetime in the 
neighbourhood of the event horizon, chosen so that surfaces of constant x° = v intersect 
the horizon on spacelike slices £„. The coordinates x M used in ( |2.2j ) provide an example of 
such a coordinate chart (as we will see these are valid over a much larger range than the 
neighbourhood of the horizon). 

As surfaces of constant v are spacelike, the null geodesies that generate the event 
horizon each intersect any of these surfaces exactly once. Consequently, we may choose 
the coordinate v as a parameter along geodesies. Then we can label the geodesies by 
a a , the value of x a at which the geodesic in question intersects the surface v — 0. The 



coordinate system {v, a a } is of the form described in § 3.1 ; as a result in these coordinates 
the entropy (d — l)-form is given by ( |3.2|) . We will now rewrite this expression in terms 
of the coordinates x^ at v = 0; for this purpose we need the formulas for the change of 



14 This definition is consistent with the Noether charge derivation of entropy currents, a la Wald, c/., 
p3| for a discussion for dynamical horizons. We review the connection with Wald's construction briefly in 
Appendix ^|. 

15 We assume here that the null energy condition is satisfied. This is true of the Lagrangian used in (33] 



to construct the gravitation background (|2.2| 
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coordinates from x M to {v, a a }, in a neighbourhood of v = 0. It is easy to verify that 



x a = a a + — v + ^— n" d„ ( — ) + 0(v d ) 
n v 2n v M \n v 1 

dx a = da a + vda k d k (—) +dv [ — + —n^d, [ — ) ) + 0(v 



(3.5) 



2x 



The coordinate transformation ( |3.5|) allows us to write an expression for the metric on 
the event horizon in terms of the coordinates {v,a a }, in a neighbourhood of v — 0. Let 
dx^ dx v = Gmn dx M dx N \n denote the metric restricted to the event horizon in the x M 
coordinates. 



dsL = H fll/ (x) dx 11 dx v = g a b da a da b 



hij ( v, a % + — J ida l + v da k d k ( — 1 J I da 7 + u da k d k I — ) ) + 0(w 



(3.6) 



where hij(v,x) is the restriction of the metric H^ u onto a spatial slice £„, which is a 
constant-?; slice. Note that since the horizon is null, all terms with explicit factors of dv 
cancel from (|3.6|) in line with the general expectations presented in § |3.1| . It follows that 
the determinant of the induced metric, y/g of fl3.2|) , is given as 

^=^h + l- Lj diVh + Vhn v d~^j + 0(v 2 ) , (3.7) 

where h is the determinant of the metric on £„, in x M coordinates (restricted to w = 0). 
We are now in a position to evaluate the area (d — l)-form 

\/h 

a = — — da 1 Ada 2 ...A da d ~ x , (3.8) 



AG 



d+l 



at v = 0. Clearly, for this purpose we can simply set to zero all terms in ( |3.5| ) with explicit 
powers of v, which implies that da a = dx a — ^ dv and 



Vh 

a 



d-l 



. dx 1 A dx 2 ... A - — dX A dx 1 A . . . A dx 1 ' 1 A dx i+1 A . . . A dx 



d-l 



(3.9) 

From (|3.9j) we can infer that the area-form can be written in terms a current as 



a = £ / M , 1M2 -^ JZ 1 dx^ A ... A daf* (3.10) 
(d-l)! 
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where Jg is given by 

_ Vh 

J s~ (d+1) nV V-^> 



and our choice of orientation leads to e„i2...(d-i) = 1- in Appendix [CJ we show that one 
can obtain this expression using the construction of an entropy (d — l)-form due to Wald, 
see ( p.12 ). We can further establish that 



da = (d - 1)! Wa - Md daJ * dx ^ A A ^ (3 ' 12) 

so that da is simply the flat space Hodge dual of Jg. While the appearance of the flat 
space Hodge dual might be puzzling at first sight, given the non-flat metric on 7i, its origins 
will become clear once we recast this discussion in terms of the fluids dynamical variables. 

3.3 Properties of the area-form and its dual current 

Having derived the expression for the area-form we pause to record some properties which 
will play a role in interpreting Jg as an entropy current in hydrodynamics. 

Non-negative divergence: Firstly, we note that the positivity of da (argued for on 
general grounds in § |3.1|) guarantees the positivity of <9 M Jg; hence we have <9 M Jg > 0. This 
in fact may be verified algebraically from (|3.7|) , as 



^—d v (^)=d^ s . (3.13) 
The positivity of d v (^fg) thus guarantees that of d^Jg as is expected on general grounds. 



Lorentz invariance: The final result for our entropy current, ( 3.12|) , is invariant under 



Lorentz transformations of the coordinate (a physical requirement of the entropy current 
for relativistic fluids) even though this is not manifest. We now show that this is indeed 
the case. 

Let us boost to coordinates = ty^x v \ denoting the horizon metric in the new coor- 
dinates by h^ v and the boosted normal vector by n M we find 

A v n m 

= A m A? h mn , A™ = A r - (3.14) 
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(where we have used h^h^ = ). It is not difficult to verify that 

n v n v 

from which it follows that ^ = thereby proving that our area- form defined on the a 
spatial section of the horizon is indeed Lorentz invariant. 



4 The Horizon to Boundary Map 

4.1 Classification of ingoing null geodesies near the boundary 

Our discussion thus far has been an analysis of the causal structure of the spacetime 



described by the metric in (|2.2|) and the construction of an area-form on spatial sections 
of the horizon in generic spacetimes. As we are interested in transporting information 
about the entropy from the horizon to the boundary (where the fluid lives), we need to 
define a map between the boundary and the horizon. The obvious choice is to map the 
point on the boundary with coordinates to the point on the horizon with coordinates 
(ruix), x^). More geometrically, this corresponds to moving along the geodesies x M = 
constant. However, congruences of null geodesies shot inwards from the boundary of AdS 
are far from unique. Hence, we digress briefly to present a characterization of the most 



general such congruence. In § |4.2| we will then see how the congruence of geodesies with 
constant x M fits into this general classification. 

We will find it simplest to use Fefferman-Graham coordinates to illustrate our point. 
Recall that any asymptotically AdS^+i spacetime may be put in the form 

2 du 2 + U]av + u d 4>^u{w)) dw^ 1 dw u 

ds = i - '- , 4.1 

u z 

in the neighbourhood of the boundary. The collection of null geodesies that intersect the 
boundary point (w^,u = 0) are given by the equations 



dw A 

~lx =u 



2 (t A + 0(u d )) (4.2) 



where A runs over the d+ 1 variables {u, u> M } and the null tangent vector must obey t A = 
0. It is always possible to re-scale the affine parameter to set t u = 1; making this choice, our 
geodesies are labelled by a d- vector t M satisfying r] fiu t^t u = —1. With these conventions 
f 1 may be regarded as a d- velocity. In summary, the set of ingoing null geodesies that 
emanate from any given boundary point are parameterized by the d — 1 directions in 
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which they can go - this parameterization is conveniently encapsulated in terms of a unit 
normalized timelike d- vector t M which may, of course, be chosen as an arbitrary function of 
x^. Consequently, congruences of ingoing null geodesies are parameterized by an arbitrary 
c?-velocity field, t^(x) on the boundary of AdS. 

4.2 Our choice of t^(x) 

It is now natural to ask what is for the congruence defined by x M = const in the 

coordinates of f33fl . The answer to this question is easy to work out, and turns out to be 
satisfyingly simple: for this choice of congruence, t^(x) = u^(x) where u^(x) is the velocity 
field of fluid dynamics^ 

While metrics dual to fluid dynamics are automatically equipped with a velocity field, it 
is in fact also possible to associate a velocity field with a much larger class of asymptotically 
AdS spacetimes. Recall that any such spacetime has a boundary stress tensor For 
most such spacetimes there is a natural velocity field associated with this stress tensor; 
the velocity u^{x) to which one has to boost in order that T 0t vanish at the point x. More 
invariantly, u^(x) is chosen to be the unique timelike eigenvector of the matrix T^(x).0 
That is, we choose u^{x) to satisfy 

(v^ + u^u u )T" K u K = (4.4) 

This definition of u^(x) coincides precisely with the velocity field in |33| (this is the so- 
called Landau frame). The null congruence given by t^(x) = u^(x) is now well defined 
for an arbitrary asymptotically AdS spacetime, and reduces to the congruence described 
earlier in this section for the metrics dual to fluid dynamics. 



Hi 



In order to sec this note that 



dx 11 dw^ du 

u„.—r— = u u 



M dX M dX dX / a o\ 

dX ~ ^ dX 



whereas indicated quantities on the LHS of (4.3) refer to the coordinate system of |33J, the quantities on 



the RHS refer to the Fcffcrman-Graham coordinates ( J4 . 1| ) . It follows from these formulae that the geodesic 



with t A = (1,1*^) maps to the null geodesic = in the coordinates used to write (2.2) 



17 In a general coordinate system the stress tensor is proportional to the extrinsic curvature of the bound- 
ary slice minus local counter-term subtractions. In the Feffcrman-Graham coordinate system described 
above, the final answer is especially simple; oc (jy^ix 11 ). 

18 This prescription breaks down when goes null - i.e., if there exist points at which the energy moves 
at the speed of light. 
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4.3 Local nature of the event horizon 



As we have seen in § ^ above, the event horizon is effectively local for the metrics dual to 
fluid dynamics such as (|2.2f) . In particular, the position of the event horizon th{x^) depends 
only on the values and derivatives of the fluid dynamical variables in a neighbourhood of x M 
and not elsewhere in spacetime. Given the generic teleological behaviour of event horizons 
(which requires knowledge of the entire future evolution of the spacetime), this feature of 
our event horizons is rather unusual. To shed light on this issue, we supply an intuitive 
explanation for this phenomenon, postponing the actual evaluation of the function rjj(x M ) 

to §0- 

The main idea behind our intuitive explanation may be stated rather simply. As we have 



explained above, the metric of [33] is tube- wise well approximated by tubes of the metric 
of a uniform black brane at constant velocity and temperature. Now consider a uniform 
black brane whose parameters are chosen as = (—1,0,0,0) and b = 1/(ttT) = 1 by a 
choice of coordinates. In this metric a radial outgoing null geodesic that starts at r = 1 + 5 
(with 5 ^> e) and v = hits the boundary at a time Sv = J pjo ~ — 4 In 5. Provided 



this radial outgoing geodesic well approximates the path of a geodesic in the metric of []33 
throughout its trajectory, it follows that the starting point of this geodesic lies outside the 
event horizon of the spacetime. 

The two conditions for the approximation described above to be valid are: 

1. That geodesic in question lies within the tube in which the metric of |33j is well 



approximated by a black brane with constant parameters throughout its trajectory. 
This is valid when 5v ~ —4 hx5 -C 1/e. 



2. That even within this tube, the small corrections to the metric of [|33| do not lead to 
large deviations in the geodesic. Recall that the radial geodesic in the metric of |33 
is given by the equation 

dv G rv +0(e) 2 + (e) 



dr G vv + 0(e) f( r ) + 0(ey 

This geodesic well approximates that of the uniform black brane provided the O (e) 
corrections above are negligible, a condition that is met provided f(r) 3> e, i.e., when 
|r - 1| = 5 > e. 

Restoring units we conclude that a point at r = r(l + 5) necessarily lies outside the 
event horizon provided 5 ^> e (this automatically ensures 8v ~ —4 ln<5 1/e. when e is 
small) . 

In a similar fashion it is easy to convince oneself that all geodesies that are emitted from 
r = |(1 — 8) hit the singularity within the regime of validity of the tube approximation 

20 



provided 5 ^> e. Such a point therefore lies inside the event horizon. It follows that the 
event horizon in the solutions of [33) is given by the hypersurface r = ttT (1 + O (e)). 



5 Specializing to Dual Fluid Dynamics 

We will now proceed to determine the precise form of the event horizon manifold to second 
order in e using the results obtained in § [| This will be useful to construct the entropy 
current in the fluid dynamics utilizing the map derived in § [|. 



5.1 The local event horizon dual to fluid dynamics 



The metric dual to fluid flows given in JHJ takes the form (|2.2|) with explicitly determined 
forms of the functions in that metric (see Appendix We list the properties and values 
of these functions that we will need below:Q 
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where C is the Catalan number. We encounter here various functions (of the boundary 
coordinates) which are essentially built out the fluid velocity and its derivatives. These 
have been abbreviated to symbols such as S3, 61, etc., and are defined ( |A.3| ). Likewise 
B°° and B fin are defined in fl02p . 

Using the equation for the conservation of stress tensor (d^T^ = 0) up to second order 
in derivatives one can simplify the expression for r H ( |2.16|) . Conservation of stress tensor 
gives 



d v 



I"' 



b 3 



(5.2) 



Projection of (~ |5.2|) into the co-moving and transverse directions, achieved by contracting 



19 Since we r equire only the values of the functions appearing in the metric (^3) and ( |2.4[ ) at r = 1/6 to 
evaluate ( [2.16 ), we present here the functions evaluated at this specific point. The full expressions can be 
found in Appendix |A]. see (A. 5) and (A. 10). 
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it with and P^ u respectively, we find 



s£ ) -^u»d,b=^a^a^ = 0(e 2 ) 

p»» ( fe 2 ^1) + = J? pu _ 3 U P df3U ^ + ( e 3) 



(5.3) 



Inserting Q5.3| ) into (2.14) we see that r^) of ( |2.14| ) simply vanishes for the spacetime dual 

- x- At next order this formula is 





to fluid dynamics, and so, to first order in e, r H 
corrected to 



,(2) 



(5.4) 



H ~ b(x) 

In order to get this result we have substituted into ( [2.16| ) the first of (p^3|), utilized the fact 
that r(i) = and the observation (from the second line of ( |5.3j )) that 



P*" {b 2 £ ] + d,b) (b 2 j^ + dv b) = 0(e 4 



In this special case the components of normal vector in the boundary directions ( [2.18 ) 
(accurate to O (e 2 )) are given by 



n 



" = (l + 4 2 ))^--P^ (d a a ai/ -3a ua u^d l3 u a )+b 2 P^j^ . (5.5) 



5.2 Entropy current for fluid dynamics 



We will now specialize the discussion of § |3.2| to the metric of |53] , using the formulae 



derived in § |5.1|. In the special case of the metric of [33] we have 



= p (l - J F 2 <V ^ + 3 ^ (2 ) + 4 2) ) 



1 

fe3 



1 



-F 2 a m a*" + U -o, v o»" + 4 2) + 4 2) , 



3_6_ 2 
T 



(5.6) 



where the various quantities are defined in (|5.1|) . We conclude from ( p.llj ) that 



r»( 5 ) 7,3 _ , u I i 6 ci2_ -.a/3, &2 _a/3 , 3 6 (2) , (2) 
4 Cj^y J g — U y 1 — r (T a/3 (T + CT + — S b + S y a ' 



+ 6 2 



l(d a a av -3a ua uPdpu a )+jW 



(5.7) 



This is the expression for the fluid dynamical entropy current which we derive from the 
gravitational dual. 
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6 Divergence of the Entropy Current 

In previous sections, we have presented a gravitational construction of an entropy current 
which, we have argued, is guaranteed to have non- negative divergence at each point. We 
have also presented an explicit construction of the entropy current to order e 2 in the 
derivative expansion. In this section we directly compute the divergence of our entropy 
current and verify its positivity. We will find it useful to first start with an abstract 
analysis of the most general Weyl invariant entropy current in fluid dynamics and compute 
its divergence, before specializing to the entropy current constructed above. 



6.1 The most general Weyl covariant entropy current and its divergence 



The entropy current in ci-dimensions has to be a Weyl covariant vector of weight d. We 
will work in four dimensions (d = 4) in this section, and so will consider currents that are 
Weyl covariant vector of weight 4. Using the equations of motion, it may be shown that 
there exists a 7 dimensional family of two derivative weight 4 Weyl covariant vectors that 
have the correct equilibrium limit for an entropy current. In the notation of (reviewed 
in Appendix [B]) , this family may be parameterized as 



(4vr r])" 1 j£ = 4 Gy b 3 Jg = [l + b 2 (A x a Q/3 a af3 + A 2 cu a p u af} + A 3 K )] 

+ b 2 [BxVxt^ + B 2 V x u» x ] 
+ Ci b t + C 2 b 2 u x V x t + ... 



a' 



(6.1) 



where b = (ttT)^ 1 ,!] = b 3 )" 1 and the rest of the notation is as in [fH| (see also 

Appendix |A| and Appendix ||). 

In Appendix [TJ we have computed the divergence of this entropy current (using the 
third order equations of motion derived and expressed in Weyl covariant language in 
Our final result is 



AGf b 3 v,r s 



b 

2 



a 



b 2Ax + 4A 3 - 



4 



In 2 u x V x a^ + 4 b (A 2 + A 3 ) u» a u a u 



+b (4 A 3 - -)(a^ a a v ) + b C 2 V»t 
(B 1 - 2A 3 ) b 2 V„V x a^ + (d + C 2 ) b 2 ^V x a» x + ... 



(6.2) 



Note that the leading order contribution to the divergence of the arbitrary entropy 
current is proportional to a^ v . This term is of second order in the derivative expansion, 
and is manifestly non-negative. In addition the divergence has several terms at third order 
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in the derivative expansion. 

Within the derivative expansion the second order piece dominates all third order terms 
whenever it is nonzero. However it is perfectly possible for to vanish at a point - a^ v 
are simply 5 of several independent Taylor coefficients in the expansion of the velocity field 
at a point (see Appendix [D] for details). When that happens the third order terms are 
the leading contributions to T>^Jg. Since such terms are cubic in derivatives they are odd 
orientation reversal (x M — > — x M ), and so can be non-negative for all velocity configurations 
only if they vanish identically. We conclude that positivity requires that the RHS of (|6.2j ) 
vanish upon setting cr M1/ to zero. 

As is apparent, all terms on the first two lines of ( |6.2|) are explicitly proportional to 
a^y. The two independent expressions on the third line of that equation are in general 
nonzero even when a^ v vanishes. As a result T>^Jg > requires that the second line of 
( |6.2|) vanish identically; hence, we obtain the following constraints on coefficients of the 
second order terms in the entropy current 

B 1 = 2A 3 C 1 + C 2 = (6.3) 

for a non- negative divergence entropy current. 

These two conditions single out a 5 dimensional submanifold of non-negative divergence 
entropy currents in the 7 dimensional space ( |3.1D of candidate Weyl covariant entropy 
currents. 

Since a local notion of entropy is an emergent thermodynamical construction (rather 
than a first principles microscopic construct), it seems reasonable that there exist some 
ambiguity in the definition of a local entropy current. We do not know, however, whether 
this physical ambiguity is large enough to account for the full 5 parameter non uniqueness 
described above, or whether a physical principle singles out a smaller sub family of this 
five dimensional space as special. Below we will see that our gravitational current - which 
is special in some respects - may be generalized to a two dimensional sub family in the 
space of positive divergence currents. 

6.2 Positivity of divergence of the gravitational entropy current 

It may be checked (see Appendix [FJ) that our entropy current (|5.7| ) may be rewritten in 



the form ( |6.1| ) with the coefficients 



1 7T In 2 „ 1 ,1 

Bl = i; B 2 = i < 6 - 4 > 
1 4' 2 

Ci = C 2 = 
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It is apparent that the coefficients listed in ( |6.4j ) obey the constraints of positivity ( |6.3| ). 



This gives a direct algebraic check of the positivity of the divergence of (|5.7|) . 

The fact that it is possible to write the current Q5.7p in the form fl6.1| ) also demonstrates 
the Weyl covariance of our current ( p.7| ). 

6.3 A two parameter class of gravitational entropy currents 

As we have seen above, there exists a five parameter set of non-negative divergence con- 
formally covariant entropy currents that have the correct equilibrium limit. An example 



of such a current was first constructed in 40 



Now let us turn to an analysis of possible generalizations of the gravitational entropy 
current presented in this paper. Our construction admits two qualitatively distinct, rea- 
sonable sounding, generalizations that we now discuss. 

Recall that we constructed our entropy (d — l)-form via the pullback of the area-form on 
the event horizon. While the area-form is a very natural object, all its physically important 
properties (most importantly the positivity of divergence) appear to be retained if we add 
to it the exterior derivative of a (d — 2)-form. This corresponds to the addition of the 
exterior derivative of a (d — 2)-form to the entropy current Jg. Imposing the additional 
requirement of Weyl invariance at the two derivative level this appears to give us the 
freedom to add a multiple of t V\uj Xct to the entropy current in four dimensions. 

In addition, we have the freedom to modify our boundary to horizon map in certain 
ways; our construction of the entropy current (|5.7|) depends on this map and we have 
made the specific choice described in § |j. Apart from geometrical naturalness and other 
aesthetic features, our choice had two important properties. First, under this map rff(x M ) 
(and hence the local entropy current) was a local function of the fluid dynamical variables 
at x^. Second, our map was Weyl covariant; in particular, the entropy current obtained 
via this map was automatically Weyl covariant. We will now parameterize all boundary 
to horizon maps (at appropriate order in the derivative expansion) that preserve these two 
desirable properties. 

Any one to one boundary to horizon map may be thought of as a boundary to boundary 
diffeomorphism compounded with the map presented in § |j. In order to preserve the lo- 
cality of the entropy current, this diffeomorphism must be small (i.e., of sub-leading order 
in the derivative expansion). At the order of interest, it turns out to be sufficient to study 
diffeomorphisms parameterized by a vector 5( that is of at most first order in the deriva- 
tive expansion. In order that our entropy current have acceptable Weyl transformation 
properties under this map, S( must be Weyl invariant. Up to terms that vanish by the 
equations of motion, this singles out a two parameter set of acceptable choices for 5(; 

= 2 8X 1 b + 5X 2 b 2 (6.5) 
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To leading order the difference between the (d — l)-forms obtained by pulling the area 
(d — l)-form a back under the two different maps is given by the Lie derivative of the 
pull-back s of a 

6s = C s< ;s = d(5( fl s^ + 5( ll (dsr. 
Taking the boundary Hodge dual of this difference we find 



5J% = £ S (Jg- JsVJC 

Similarly 



v v [jg 5C - J u s <*ci + <*C v v r s (6 ' 6) 



Using the fluid equations of motion it turns that the RHS of ( |6.6| ) is of order e 3 (and so zero 
to the order retained in this paper) for oc b 2 Z M . Consequently, to second order we find a 
one parameter generalization of the entropy current - resulting from the diffeomorphisms 
( |6.5D with 6X2 set to zero. 

Note that, apart from the diffeomorphism shift, the local rate of entropy production 
changes in magnitude (but not in sign) under redefinition ( |6.6| ) by a factor proportional 
to the Jacobian of the coordinate transformation parameterized by 5(. In Appendix [B] wc 
have explicitly computed the shift in the current (|5.7|) under the operation described in 
( |6.6D (with 5( of the form ( |6.5|) ) and also explicitly verified the invariance of the positivity 
of divergence under this map. 

In summary we have constructed a two parameter generalization of our gravitational 
entropy current (|5.7|) . One of these two parameters arose from the freedom to add an exact 
form to the area form on the horizon. The second parameter had its origin in the freedom 
to generalize the boundary to horizon map. 



7 Discussion 

We have demonstrated that any singularities in the metrics of |33|], dual to fluid dynamics, 
are shielded behind a regular event horizon (we expect the same to be true for the solution 
of [0). Further, we have shown that the structure of this event horizon is determined 
locally by the variables of fluid dynamics, and have presented an explicit expression for 
the location of the event horizon to second order in the e (boundary derivative) expansion. 
Remarkably, the event horizon, which is a global concept in general relativity, turned out to 
be rather simple to locate, partly due to our choice of particularly useful coordinate system 
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(O), and more importantly due to the long- wavelength requirement that our solution be 
dual to a system described by fluid dynamics. We emphasize that within the boundary 
derivative expansion of this paper we are directly able to construct the event horizon; 
we did not need to discuss other more local constructs like the apparent horizon as an 
intermediate step towards understanding the global structure of our solutions. 

We have also constructed an entropy (d — l)-form on the event horizon of an arbitrary 
d+1 dimensional spacetime and used the pullback of this form to the boundary to construct 
a manifestly non-negative divergence entropy current for asymptotically AdS^+i solutions 
of gravity with a horizon. We have derived an explicit expression for this entropy current 
for the solutions dual to and demonstrated a direct algebraic check of the positivity of 
divergence of this current within fluid dynamics. 

In order to lift the entropy (d — l)-form from the horizon to the boundary, we used a 
natural map between the horizon and the boundary, given by ingoing null geodesies which 
emanate from the boundary in the direction of the fluid flow. These ingoing geodesies in 
fact determine the coordinate system of |3l| (they constitute lines of constant x M for the 
metric Q). 

We also directly studied a seven parameter family of weight four Weyl covariant fluid 
dynamical vectors that have the appropriate equilibrium limit to be an entropy current and 
demonstrated that a 5 parameter subclass of this family of currents has non negative di- 
vergence to second order in the derivative expansion. The entropy currents we constructed 
via a pullback of the area form constitute a special subclass of these currents. It is natural 
to inquire what the gravitational interpretation of the remaining currents is.0 It is natural 
to wonder whether they are assosicated with appararent horizons^] and other quasi-local 
horizons (such as trapping/dynamical horizons, isolated horizons, e.g.,[[|4], ^5[). At least 



several of these horizons also appear to obey versions of the area increase theorem. Con- 
sequently, it should be possible to obtain conserved entropy currents via the pullback of 
a suitably defined area form on these horizons. Apparent and other dynamical horizons 
have one initially unpalatable feature; their structure depends on a choice of the slicing 
of spacetime into spacelike surfaces. However perhaps it is precisely this ambiguity that 
allows these constructions to cover the full 5 parameter set of non negative entropy currents 
discussed above ?Q Note that in the context of dynamical horizons, obtain^] a char- 
acterization of a membrane fluid obeying non-relativistic hydrodynamics equations with 
a uniquely specified entropy. Their system has rather different characteristics (absence of 



20 We thank M. Van Raamsdonk for raising this question. 

21 We thank A. Stromingcr for stressing the physical relevance of apparent horizons to our situation, and 
for a very useful related discussion. 

22 A cautionary note is in order; it is possible that a subclass of the 5 parameter non negative divergence 
entropy currents is an artifact of the derivative expansion, and has no continuation to finite e. We thank 
M. Van Raamsdonk for discussions on this point. 

23 We thank I. Booth for bringing this reference to our attention. 
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shear- viscositly for instance) and appears to model the black hole as a fluid, rather than 
construct an explicit dual as in the current discussion. It would interesting to understand 
this connection better. 

We re-emphasize that our results demonstrate that each of the solutions of (with 
regular fluid data) has its singularities hidden from the boundary by a regular event horizon. 
Consequently all gravitational solutions dual to regular solutions of fluid dynamics obey 
the cosmic censorship conjecture. It would be interesting to investigate how our results 
generalize to irregular (e.g., turbulent) solutions of fluid dynamics, as also to gravitational 
solutions beyond the long wavelength expansion. As we have explained in the Introduction, 
several such solutions are dual descriptions of the field theoretic approach towards local 
equilibrium. The appearance of a naked singularity in this approach would appear to imply 
singularities of real time correlation functions in this process. It would be fascinating to 
study this connection in more detail. On a more speculative, or perhaps more ambitious 
note, it is natural to inquire what (if any) feature of field theoretic correlators would be 
sensitive to the apparently crazy nature of near singularity dynamics even when the latter 
is cloaked by a horizon. 

Another interesting direction concerns a' corrections to the bulk solutions, which in the 
language of the dual field theory correspond to finite 't Hooft coupling effects. In the bulk, 
there is a well developed formalism due to Wald [43, 4?| which provides a generalization 
of the Bekenstein-Hawking area formula for the entropy of the black hole to higher deriva- 
tive gravity. The main idea is to construct the entropy of (asymptotically flat) solutions 
using a variational principle of the Lagrangian; essentially, from the variational principle 
one obtains the first law of thermodynamics, which is used to construct the entropy as a 
Noether charge. This construction of the Noether charge entropy is conceptually similar 
to the area-form we present and in fact reduces to it in the two derivative limit. It might 
be possible - and would be very interesting - to generalize the discussion presented in 
this paper to be able to account for a' corrections (see Appendix [Cj). The key issue here 
would be to find an analogue of the area increase theorem for a' corrected gravity. This 
is complicated from a pure general relativity standpoint, owing to the fact that higher 
derivative theories generically violate the energy conditions.^ The AdS/CFT correspon- 
dence seems to require that such a generalization exist, and it would be very interesting to 
determine it. It is possible that the requirement of the existence of such a theorem provides 
'thermodynamical' constraints to a' corrections of the low energy equations of gravity.0 



24 In the context of supersymmetric solutions in a' corrected gravity, as discussed in |4q] for the so called 
small black holes, while the entropy remains proportional to the area and specific classes of solutions satisfy 
the averaged null energy condition [fl9|| ; one still is unable to show the desired monotonicity property of 
entropy. 

25 In |Q the authors demonstrate the second law for the Einstein Hilbert action deformed by an R 2 
term; this Lagrangian however is not likely to arise as the low energy effective action from string theory 
fell. See also |52|, p3[ for recent discussions of constraints on parameters appearing in higher derivative 
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If the gravity /fluid dynamics correspondence could be understood in more detail for 
confining gauge theories (see e.g., |54], |55]|), fluid dynamics could give us a handle on very 
interesting horizon dynamics. For instance, one might hope to explore the possibility of 
topological transitions of the event horizon.^ 

Turning to more straightforward issues, it would be interesting to generalize the dis- 
cussion of this paper to encompass the study of field theory arbitrary curved manifolds. It 
would also be interesting to generalize our analysis to the bulk dual of charged fluid flows, 
and especially to the flows of extremal charged fluids. This could permit more direct con- 
tact with the entropy functional formalism for extremal black holes. A natural framework 
for such analysis, specifically in relation to the horizon dynamics studied here, is provided 



by the near-horizon metrics for degenerate horizons discussed in ||56|| . Furthermore, one 
might hope that such a study could have bearing on studying the behaviour of superfluids. 

Finally, note that while field theoretic conserved currents are most naturally evaluated 
at the boundary of AdS, the entropy current most naturally lives on the horizon. This is 
probably related to the fact that while field theoretic conserved currents are microscopically 
defined, the notion of a local entropy is an emergent long distance concept, and so naturally 
lives in the deep IR region of geometry, which, by the UV/IR map, is precisely the event 
horizon. Correspondingly, we find it fascinating that, in the limits studied in this paper, the 
shape of the event horizon is a local reflection of fluid variables. This result is reminiscent 
of the membrane paradigm of black hole physics. It would be fascinating to flesh out this 
observation, and perhaps to generalize it. 
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theories in connection to the dual hydrodynamic description. 

26 Whilc Cosmic Censorship precludes splitting of black holes, they can easily merge without the curva- 
tures becoming large. 
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A Notation 



We work in the mostly positive, ( — h + . . .), signature. The dimensions of the spacetime 
in which the conformal fluid lives is denoted by d. In the context of AdS/CFT, the dual 
AdS^+i space has d + 1 spacetime dimensions. The event horizon is a d-dimensional null 
manifold Ti. Ti is foliated by d — 1 dimensional constant v spatial slices denoted by E„. 
The induced metric on E„ is denoted by h a b (and h denotes its determinant). 

Latin alphabets A, B, . . . are used to denote the d + 1 dimensional bulk indices which 
range over {r, 0, 1, ... , d— 1}. Lower Greek letters fi, u, . . . indices range over {0, 1, . . . , d— 1} 
and lower case Latin letters a,b, . . . indices range over {1, . . . , d — 1}. The co-ordinates in 
the bulk are denoted by X A which is often split into a radial co-ordinate r and x M . We 
will often split x M into v and x a . 

In these co-ordinates, the equation for the horizon takes the form S-^ = r — r^(x) = 0. 
We can choose to eliminate the co-ordinate r in favour of x M 's via this equation. Then, 
in the co-ordinates the components of the metric are denoted by H^ v . In addition, 
we find it convenient to use a co-ordinate system a a , A on 7i - in these co-ordinates, the 
components of the induced metric on the horizon take the special form g\\ = g a \ = and 



In Table [T], we list the physical meaning and the definitions of various quantities used 
in the text, referring to the equations defining them where appropriate: 

A.l Fluid dynamical parameters 

Various expressions in the text and are built out of the fluid velocity; we list them here for 
convenience. The basic building blocks are the derivatives of the fluid velocity, decomposed 
into appropriate representations based on their symmetries. We have (see the table above 
for the physical meaning of these parameters), 



9ab ^ 0. 

Our convention for the Riemann curvature tensor is fixed by the relation 



(A.l) 



V 



4iu 




(A.2) 
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Symbol 


Definition 


Symbol 


Definition 


A 
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UlillcilolUilo Ul UUU.IlU.cLlj 


n 
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X A 


Bulk co-ordinates 


x v 


Boundary co-ordinates 




A spatial slice of Ji 




Co-ordinates on Ti 


Gab 


Bulk metric, (gg) 




Boundary metric (Minkowski) 


h a b 


Induced metric on Yj v 


9ab 


Metric on C Ti 


r H (x) 


Horizon function, (|2.7|) 




Induced metric on H 


S n = 


Eqn. of Horizon 


s 


Entropy (d-l)-form on Sa 


Z A 


Normal vector to the Horizon (|2.9|) 




See (|2.170 
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See (£J), (|A.5|), (|5J) 
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See ( 
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Jn 


See (O), (IA.10I), (O) 
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Shear viscosity 
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Entropy current 
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Projection tensor, t]^ + u^u u 


a 11 


Fluid acceleration, (|A.2) 




Fluid expansion, (A. 2) 




Shear strain rate, (|A.2|) 




Fluid vorticity, (A. 2) 




Visco-elastic stress 












Weyl-covariant derivative 


A/A 


See ( 


B.4 


) 


T3 a 


Riemann tensor 


J fJ,V 


V^A 
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V U A^ 




Ricci tensor /scalar 




See ( 


B.4 
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Einstein tensor 




See ( 


B.4 
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C 1 'fiuXcr 


Weyl curvature 











Table 1: Conventions used in the text 



In addition, we will have occasion at various points in the text to encounter various 
functions built out of the first derivatives of the fluid velocity defined in (|A.2|) . These 
functions were defined in |33| to present the second order metric, and show up for example 
in (I5TTI). We have: 



s 3 = T p«f> d a dab 6i = Vu a Vu a , 
o 
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(A.3) 



(IV 



where T> = <9 M (Note that this is a different derivative from the Weyl covariant derivative 
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introduced in Appendix [B|; the distinction should be clear from the context). 
9 



V4i 



1 -p: d, (d pUa + d aU/3 ) - ~ p ap pi d,d aU p 



-p»p^d a d^ 



V5 v = P?P af> d a d f! u ll 

23 1„ = d a u a Vu v , 2J2j, = e a f3 lv u a VvP P , 233„ = a au Vu c 



(A.4) 



A. 2 The functions appearing in the second order metric 



The metric fl2.2p derived in [33] has been rewritten in terms of various auxiliary functions 
used to define S(r, x M ) and Xfj,u(r,x^). These functions can be read off from Eq (5.25) of 
; we list them here for convenience .p 7 ] 



Scalars under 50(3) spatial rotations: The scalar functions appearing at first and 
second order are respectively, 



4 2) (r,0 = ^ 2 h 2 (b r) 



,(i) 



(r, x M ) = - r d\U 



1 k 2 (br) 

^2 



in terms of several functions of r which are given as 



F(r) 



In 



(l + r) 2 (l + r 2 ) 



2 arctan(r) + n 



Defining 

S h (r) = 



3^6 4 + iWfc(r)6 5 



(A.5) 



(A.6) 



5 fc (r) = 12r 3 fc 2 (r) + (3 r 4 - 1) ^(r) - y s 3 + 2 r 6i - ^ 63 + 64 + 5 W*(r) 65 

(A.7) 



27 One notational change wc have made is to rename the functions 01$ appearing in [^3) to tffl. We 
don't list this here as it doesn't appear directly in our analysis of the entropy current. 
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where the functions Whir) and Wk(r) are given by 



4 (r 2 + r + l) 2 - 2 (3 r 2 + 2r + 1) F(r) 



W h (r) = - 
W k (r) 



3 r (r + iy (r 2 + iy 

2 4 (r 2 + r + 1) (3 r 4 - 1) F(r) - (2r 5 + 2r 4 + 2r 3 - r - 1) 

3 r (r + 1) (r 2 + l) ' 

The other symbols S3, ©1, etc., are defined in ( |A.3|) . We can now write the expressions 

(2) 

for the functions appearing in the definition of s b as 



1 f°° dr f°° f 1 

Hr) = -^Sfr + l ^ jf dyy^S h (y)--S h ° 



(A.{ 



where we have defined 

5^=^64+^65) , fi fc °°= (-^3 + 261-^63-^64+^65) • (A.9) 

Vectors under 50(3) spatial rotations: The vector functions appearing at first and 
second order are respectively, 

1 / r 2 f°° f°° ( 1 \\ (A 10) 



where 



with 



RM _ (2r 3 + 2r 2 + 2r-3) B°° + B fi ° 

" 18r3fr + l) fr 2 + l) (A ' llj 



B°° = 4 (10 V4 + v 5 + 3 QJi - 3 232 - 6 5J3) 
B fin = 9 (20 v 4 - 5 232-6 033), 

The symbols vk and 03/c are defined above in ( |A.4|) as derivatives of the fluid velocity. 



B Weyl covariant formalism 

In this appendix, we present the various results related to Weyl covariance in hydrodynam- 
ics that are relevant to this paper. The conformal nature of the boundary fluid dynamics 
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strongly constrains the form of the stress tensor and the entropy current [f3^, [1(J. An effi- 
cient way of exploiting this symmetry is to employ a manifestly Weyl-covariant formalism 
for hydrodynamics that was introduced in the reference ||D| . 

In brief, for an arbitrary tensor with weight w, one defines a Weyl-covariant derivative 



V x ()>;;;■ = V A Q'l: 



w A x QZz 
5%A V - 6? Ax] 



(b.i; 



where the Weyl-connection A„ is related to the fluid velocity via the relation 



A^i = u x V x u^, 



V x u x 



:b.2] 



We shall exploit the manifest Weyl covariance of this formalism to establish certain results 
concerning the entropy current that are relevant to the discussion in the main text. 

B.l|, we write down the most general Weyl-covariant entropy current and compute 



In 



its divergence. This computation leads us directly to an analysis of the constraints on the 
entropy current imposed by the second law of thermodynamics. This analysis generalizes 
and completes the analysis in ||(| where a particular example of an entropy current which 
satisfies the second law was presented. Following that, in § B.2| , we rewrite the results of this 
paper in a Weyl-covariant form and show that the expression for the entropy current derived 
in this paper satisfies the constraint derived in § [B.l| . This is followed by a discussion in 
§ [B.3| on the ambiguities in the definition of the entropy current. 



B.I Constraints on the entropy current: Weyl covariance and the second law 

We begin by writing down the most general derivative expansion of the entropy current 
in terms Weyl-covariant vectors of weight 4.0 After taking into account the equations of 
motion and various other identities, the most general entropy current consistent with Weyl 
covariance can be written as: 

(4tt V y l J£ = 4 Gf b 3 J£ = [1 + b 2 (A t a aP a aP + A 2 u af3 + A 3 K ) ] u» 

+ b 2 [Bt T? A cx' lA + B 2 V x u» x ] (B.3) 
+ Cib£» + C 2 b 2 u x V x t + ... 

28 In contrast to the analysis in the main text, we find it convenient here to work with an arbitrary 
background metric, whose associated torsion- free connection is used to define the covariant derivative V^. 
29 We will restrict attention to fluid dynamics in 3 + 1 dimensions. 
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where b = (irT) 1 and we have already assumed the leading order result for the entropy 
density s = 4tt rj = (AG^ 6 3 )" 1 and £" = e a ^uj aP u v ^\ 

Now, we want to derive the constraints imposed by the second law on the A,B and C 
coefficients appearing above. To this end, we take the divergence of the entropy current 
above to get 



4 GP b 3 V^ = - 3b- 1 vP V^b -2C X FV^b 

+ b 2 V^[(A 1 a a pa af3 + A 2 u a pu af3 + A 3 'R) u» (B.5) 
+ (B 1 V x a" x + B 2 V x oo» x + C 2 u x V X F)] + ... 

where we have used the facts that T)^ = and that T>pb gets non-zero contributions only 
at second order (P-7| ). Further, u x gets non-zero contributions only at third order (the 
equations of motion force u^T^x = at second order). 

In order to simplify the expression further, we need the equations of motion. Let us 
write the stress tensor in the form 

= (16tt G^ & 4 )" 1 (jT + 4 u*V) + tt^ (B.6) 
where ir^ is transverse - u u ituv — 0. This would imply 



,3 b \ (B.7) 

4 \-urph-— a, lv 7r^ 1=0 



= b A Ufl V V T>» = b A V v {u^) - b A {V u uJ T» v 

-^iib (T„v 

b M 4t] m 

(5) 

where we have multiplied the equation by 16irG N in the second line to express things 
compactly. Similarly, we can write 2 V^b = — b 2 i^Vxcr^ which is exact upto third order 
in the derivative expansion. Note that these are just the Weyl-covariant forms of the 
equations that we have already encountered in ( [5.31 ). 

We further invoke the following identities(which follow from the identities proved in the 



30 Wc shall follow the notations of [Q in the rest of this appendix. In particular, we recall the following 
definitions 



11 = R- &V X A X + 6A X A X ; V^u v = + lo^ ( R4 ) 



Note that in a flat spacetime, R is zero but 1Z is not. Though we will always be working in flat spacetime, 
we will keep the i?-terms around to make our expressions manifestly Weyl-covariant. 
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Appendix A of |3) Q 



V^uo aP u af3 u») = AO** 00^ u au - 2 Vp x ^ x 

V^TZ m") = -2 TV" + V„ [-2 V x a» x + 2 2W A + 4 w A j^ A ] 
-2 a^-RT = 4 <v [u a D a gt^ + lu^ u a v + a a u - u a u p ] 
V,(u x V x l») = V,(£ x V x u») - tu v 
V^i x V x u») = <v VH V + £ M D x u^ x 



(B.8) 



to finally obtain 



4Gf b 3 V^ 



b 2 a 



b 2 a 



I"' 



7T 



4r/6 



2 Ai w A V x a" v + AA 2 u a u -2A 3 TV" + C 2 



(B 1 - 2A 3 ) b 2 V^V^ + (d + d) 6 2 £ M V x a 



/(A 



7T 



Arjb 



{2A l + 4A 3 ) n A P A a^ + 4 (A 2 + A 3 ) w a " + 4 A 3 a 



+ {B 1 - 2A 3 ) 6 2 T) x u^ x + (d + d) & 2 IP** + ■■■ 



d 

(B.9) 



Substituting the value of as calculated from the known stress tensor, we find 



b 2 a 



it i' 



a 



1 1 



— + I2A + 4A 3 -- + - In 2 ) u A V x a^ 



+ 4(A 2 + A 3 ) u» a + (4 A 3 - -) (a^ a a a v ) + C 2 VH V 
+ (B 1 - 2A 3 ) b 2 V,V x a» x + (d + d) &%X>a^ + . . . (B.10) 

This expression can in turn be rewritten in a more useful form by isolating the terms that 
are manifestly non-negative: 



4Gg ) & 3 2V£ 



b 
2 



2Ai + 4A 3 - 1 + i In 2 



n 2 



+ 6 (4 A 3 - -) (<r" Q t7 a ") + b d W 
+ (5! - 2A 3 ) 6 2 X? M X? A 6^ A + (d + d) & 2 



(B.11) 



The second law requires that the right hand side of the above equation be positive 



31 Since we are only interested in the case where boundary is conformally flat, we will consistently neglect 
terms proportional to the Weyl curvature in the following. 
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semi-definite at every point in the boundary. First, we note from ( |B.11| ) that the first two 



lines are positive semi-definite whereas the terms in the third line are not - given a velocity 
configuration in which the third line evaluates to a particular value, as argued in the main 
text, we can always construct another configuration to get a contribution with opposite 
sign. Consider, in particular, points in the boundary where = - at such points, the 
contribution of the first two lines become subdominant in the derivative expansion to the 
contribution from the third line. The entropy production at these points can be positive 
semi-definite only if the combination the coefficients appearing in the third line vanish 
identically. 

Hence, we conclude that the second law gives us two constraints relating A,B and C, 

viz., 

B l = 2A 3 d + C 2 = (B.12) 

Any entropy current which satisfies the above relations constitutes a satisfactory proposal 
for the entropy current from the viewpoint of the second law. 

One simple expression for such an entropy current which satisfies the above require- 
ments was proposed in |4TJ. The J x proposed there is given by 



(4vr 7])- 1 J x s = u x - h ^ [(In 2 a^a^ + cu^ cu^) u x + 2u^ (G» x + F*) + 6 V v u Xv ] 



Now, using the identity 



(B.13) 



(G^ + F X ) = -j u x - V v a Xv - V v u Xv + 2 UyF** (B.14) 

and the equations of motion, we can rewrite the above expression in the form appearing 
in (|B.3|) to get the value of A,B and C coefficients as 



A x = ; A 2 = — ; A 3 = - 

8 8' 8 

B-\ = — : Bo = — 
4' 2 

Cx = C 2 = 



(B.15) 



It can easily be checked that these values satisfy the constraints listed in ( |B.12j ). Further, 
for these values, the divergence of the entropy current simplifies considerably and we get 

4Gf 6 3 2V£ = ^a^ (B.16) 
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However, as the analysis in this section shows, this proposal is just one entropy current 
among a class of entropy currents that satisfy the second law. This is not surprising, since 
(as was noted in [HU|) the second law alone cannot determine the entropy current uniquely. 



B.2 Entropy current and entropy production from gravity 

We now calculate the coefficients A»'s and B^s for the actual entropy current calculated 
from gravity in (|5.7f) and check whether the they obey the constraints in (|B.12j ). Unlike 
the proposal in ||(| , the entropy current derived in § ^| takes into account the detailed 
microscopic dynamics(of which hydrodynamics is an effective description) encoded in the 
dual gravitational description. 

In order to cast the entropy current in the form given by ( |B.3|) , we have to first rewrite 
the quantities appearing in this paper in a Weyl-covariant form. We have the following 
relations in the flat spacetime which identify the Weyl-covariant forms appearing in the 
second-order metric of |33| - 



©4 = 2 u af3 uj af3 ; ©5 = a a p a af3 ; 
- ^ s 3 + 2 6 1 - ^ 63 = ^a a p a aP - | LU a ptu^ + l -TZ 



5 5 5 _ 5 m 11 \ (B.17) 

g v 4m + - v 5 , + - - - Qj2 M - - 5j 3 , = p; zW 

15 1 „ 1 „ 1 



y v 4m - - v 5m - - + - 5J 3 m = P v , V x uJ„ x 

These can be used to obtain 

b^ = i8p;d a ^ + i8p;p a ^ a 



A 



= 18 {-Gap a aP + to a p cu al3 ) m m + 18 V x a x + 18 V x u„ 
Bj n = 54P;V x a x u +90P»V x uj x u 

= (-54 a a p a a<i + 90 u a p uj ap ) + 54 £> A <r/ + 90 V x uj 



;b.is) 



A 



Hence, all the second-order scalar and the vector contributions to the metric can be written 
in terms of three Weyl-covariant scalars <J a p<J a P, u a pu a ^ and 1Z and two Weyl-covariant 
vectors T) x o^ and T> x uj^ . 

Using the above expressions, we can rewrite the second order scalar and the vector 
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contributions to the entropy current appearing in ( |5.1| ) as 



,(2) 



,(2) 



,(2) 



1 2 



1 57T 2 

+ ln2 + 3C + ? + - 



3, 7T 

2"' 2+ 4 



a(3 



IT 5 7T" 

_ + _ hl2 + c + _ + _ 



a a/3 a 



af3 



■ uj a/3 u af3 + - K 



1 

6 



(B.19) 



while the vector contribution is given as 



J ft 



(2) _ pv 



- V x o x + - V x u x 



] u^-V x o x + -V x uj x 

(B.20) 

Now, we use Q5.4|) , ( |5~5"1) and ( |5.6|) to write rff,n^ and in Weyl covariant form as 
follows: 



r » = b 1 + T 



5 2, , 7T 5r 
- + -n2 + C + - + 

6 3 6 48 



a/3 



2 6 



(B.21) 



b 2 



1 ^ ?r 5ir 2 

_ +ln2+3C+i + _ 



3 , 7T 

- In 2 + - 

2 4 



a/3 



1 



+ b 2 p; 



- V x a x + - V x u u x 



(B.22) 



r- 1 (, b " 
1+ 4 



(2 + In 2 + ~) a a(3 a a ? - ^ cu aP u^+^TZ )j (B.23) 
Putting all of these together we can finally obtain the expression for the entropy current: 



4Gf 6 3 J£ 









7T 




[(§- 


- in 2 H 








4 


16. 



CT — LO a pUJ ' + - It 



U> 



2 p; [\v x a x + l -v x u u x 



>a(3 1 



9 (^V x a» x + 1 -V x uj^ 







1 1 7T \ 


(l + b 2 


ig- 


- - In 2 + — 




4 16/ 



'a/3 1 



(B.24) 



from which we can read off the coefficients A, B and C appearing in the general current 
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(ED 



1 7T In 2 

Al = 4 + i6 + ^r ; 



A 9 



-4, 



B-\ = — ; Bo = — 
1 4' 2 

Ci = C 2 = 



(B.25) 



These coefficients manifestly obey the constraints laid down in (|B.12j ) and hence, the 
entropy current derived from gravity obeys the second law. Further, we get the divergence 
of the entropy current as 



AG { $ b 3 J» = b 2 a, u 



a 



1 7T 



— + 2 (- + — + - In 2) w A £> A cx^ 
2b V 4 16 



+ . . . 



b 
2 



( ~ + ^ + | In 2 



(B.26) 



which can alternatively be written in the form 
T2Vs = 2i, 



^ +4+61 " 2 W<- 



+ . . . 



(B.27) 



which gives the final expression for the rate of entropy production computed via holography. 



B.3 Ambiguity in the holographic entropy current 

We now examine briefly the change in the coefficients A, B and C parametrizing the 
arbitrary entropy current, under the ambiguity shift discussed in § see Eq. (JK.fip. In 
particular, we want to verify explicitly that under such a shift, the entropy production still 
remains positive semi-definite. 

The first kind of ambiguity in the entropy current arises due to the addition of an exact 
form to the entropy current. The only Weyl covariant exact form that can appear in the 
entropy current at this order is given by 

4 G% ] b 3 5J% = 5X b 2 V u uj^ (B.28) 



which induces a shift in the above coefficients B 2 — > B 2 + SXo. 

The second kind shift in the entropy current (due to the arbitrariness in the boundary 
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to horizon map) is parametrised by a vector (which is Weyl- invariant) and is given by 



S C 5 5 (B.29) 

= v„[j»6e-r s 6C] + 6ev u r s 

where in the last line we have rewritten the shift in a manifestly Weyl-covariant form. 
If we now write down a general derivative expansion for ££ M as 

= 2 5X l bu fl + 5X 2 b 2 F + . . . (B.30) 

the shift in the entropy current can be calculated using the above identities as 

4 G { p b 3 6Jg = 5\! b 2 a a p a a(3 u> 1 + ... (B.31) 

which implies a shift in the above coefficients given by A\ — > A\ + b\\ . 

Note that both these shifts maintain the constraints listed in (B.12| ) and hence, the 
positive semi-definite nature of the entropy production is unaffected by these ambiguities 
as advertised. 



C Wald's "entropy form" 

In this section we briefly discuss the notion of a local "entropy form" , as defined by Wald, 



47] , |57| . This is defined using a variational principle for any diffeomorphism invariant 
Lagrangian £ to derive an expression for the first law of black hole mechanics. We consider 
a d + 1 dimensional spacetime with metric Gab which is a solution to £'s equations of 
motion and denote to be the associated covariant derivative. 

C.l Stationary black branes 

Let us first consider the case of a stationary black brane, characterized by a Killing horizon 
7i which is generated by a Killing vector \ A - We normalize \ A by the condition that it 
satisfies xa V a x B — X B 011 W and assume that 7i possesses a bifurcation surface E&. 
Consider the following d — 1-form 



2i^J^G dC 



It has been shown in H3I, O that the entropy of the black hole S is then simply the integral 



of S over £&, and it satisfies the first law of thermodynamics. Hence, flC.lf) provides a local 
expression for the entropy- form. As is to be expected, this expression is not unique, and 
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suffers from ambiguities arising from: (i) the possibility of adding exact derivatives to £, 
(ii) addition of a (d — 1) form to S which arises from the additive ambiguity of the Noether 
current up to the Hodge dual of an exact d-form, and (hi) the possibility of adding to S 
an exact d — 1 form without changing the entropy S, cf. proposition 4.1 of [I3|. However, 
in the discussion that follows, these additional terms will not be important. 

It is easy to evaluate the above expression ( p.lj ) in case of General Relativity. In this 
case £ = - \ d+1) (R + A), and 

ID7T G » 



N 



or f £ = ^— (G AC G BD - G BC G AD ). (C.2) 

dR AB cD 32ttG^ +1) 

Further, on £ b , 

V[aXb\ = kab (C.3) 

where yiab is the binormal to £&, defined by 

n A B = N aXb - N bXa , (C.4) 

where N A is the "ingoing" future-directed null vector, normalized such that N A xa — — 1 
0. It is easy to show that the volume element on E& is given by (see Eq. (12.5.34) of |58j): 



= ~ 2{d--l)\ eAl - A ^ AB ~ ~ (d--l)\ €m - a <i-^ ab V cXd G ac G bd (C.5) 
Putting all this together, the entropy form S becomes 

Sal..^ = (d+1) ^...Aw (C.6) 

4C, JV 

The total entropy is given by the integral 



S = I S^.A^cfar 41 A ... A dx Ad ^ = ^ {d+l] Area (E 6 ) (C.7) 



AG N 



thus reproducing the usual Bekenstein-Hawking formula. For stationary black holes the 
area of the bifurcation surface of course coincides with the area of the black hole horizon. 

To be explicit, let us consider the example of the five dimensional stationary (boosted) 
black brane solution which is given by Q2.2p , with e = and b and m m constants (independent 



32 This normalization, together with fact that N A , \ A are both normal to £& uniquely fixes N A . 
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of The horizon is located at r = ru = 1/6- Let us consider a spacelike slice £ C Ti 
defined by dx^ = 0. The binormal uab to this surface (|C.4| ) is given in terms of the 
null vectors;^ and N A . We have the normalized Killing vector X A q^x = ^''^i and 
N A dX A = K(2dr-r 2 f(br) u^dnf). Here k = \ (r 2 d r f(br)) \ r=rH . 

From ( |C4| ) we find that the only non- vanishing components of the entropy (d — l)-form 
S are given by 



„<i-l 



We have used the fact that on the (d — l)-surface E, \fh = r d H 1 u v . 

The entropy form a, given in ( |3.10| ), agrees with the above expression: 



a = S ^dx* 11 A dx» 2 A ... A dx^- 1 (C.9) 



where we note that the vector n M in ( 3-10 ) becomes equal to in the static case (we will 



discuss the dynamical situation below). 

It is interesting to note that in case of higher derivative gravity, the entropy form 
has terms in addition to the area-form. For example, in case of Lovelock gravity, with 
Lagrangian density 

£ = R + a ( r abcd R ABCD - 4 Rab R AB + R 2 ) , (CIO) 

the entropy form S as defined by ( |C.1|) (see Eq. (72) of |[43|| ) leads to an entropy, which 
has contributions from the first Chern class of the bifurcation surface: 

S= L_ y Area(S 6 ) +8vra ! R {d ~ 1] y/^V) d d ~ l x. (C.ll) 

4:G N Js b 

C.2 Dynamical horizons 

The horizon of dynamical black holes (such as the generic situation with ( |2.2| )) is not 
generated by a Killing field and generically one doesn't have a bifurcation surface. So the 
formula ( |C.1|) cannot be applied as such. However, as argued in j43], the simplest way 



to proceed in this case is to develop a notion of a local bifurcation surface, and construct 
a "local Killing field" %■ m case °f General Relativity, this leads to a definition of the 
entropy {d — 1) form S as in ( |C.1|) ; the main distinction is that VcXd is interpreted as 
the binormal iicd only in a sufficiently small neighbourhood close to the initially chosen 
surface. Hence locally, we can continue as before using the result (|C.5|) to arrive at the 
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expression 



4G 



(d+i) 

N 



(C.12) 



Here in constructing the binormal ( |C.4| ) we have used the fact that \ A oc n A on where 



n 



A ■ 



is defined by ( |2.17| ). As explained in [|3], at any point p on a spacelike surface a C H 



we can choose coordinates such that the expression ( |U.1| ) for the entropy remains correct 
(in particular, one can choose to ensure that the additional terms arising from ambiguities 
in defining S vanish), so that the above derivation of ( |(J.12[ ) remains valid. It is easy to see 
that the entropy form a, given in ( ^.10[ ), agrees with the above expression in the dynamical 
case as well. 



C.3 Second law 

We saw above, in case of General Relativity in arbitrary dimensions, that the Wald def- 
inition of entropy leads to the area-form on the horizon. The divergence of the entropy 
current therefore is nonnegative as a consequence of Hawking's area theorem and hence 
obeys the second law of thermodynamics (assuming cosmic censorship). Recall that area 
theorem requires that the energy conditions hold; physically, the only when gravity is at- 
tractive are we guaranteed area increase. However, in case of higher derivative theories, it 
is not clear whether the second law is obeyed |I3|, [47] by the Wald entropy (c/., |50|, |59| 



for a discussion in certain special classes of higher derivative theories). This is simply be- 
cause higher derivative theories violate the energy conditions and the situation is further 
complicated by the fact that the entropy starts to depend on the intrinsic geometry of the 
black hole horizon. On the other hand, from the viewpoint of the boundary theory, a' 
corrections simply provide a one-parameter deformation of various parameters of the fluid 
which must continue to obey the second law of thermodynamics. It would be interesting 
to resolve this puzzle (see § [7| for comments). 



D Independent data in fields up to third order 

There are 16, 40 and 80 independent components at first, second and third orders in 
the Taylor expansion of velocity and temperature. These pieces of data are not all 
independent; they are constrained by equations of motion. The relevant equations of 
motion are the conservation of the stress tensor and its first and second derivatives^ (at 

33 For each independent function we count the number of independent partial derivatives at a given order; 
for the temperature we have <9pT, d^d^T, etc.. 

34 The relevant equations are just the moments of the conservation equation which arise as local con- 
straints at higher orders. 
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our spacetime point) which are 4, 16 and 40 respectively in number.^ The terms that 
appear in the three kinds of equations listed above start at first, second and third order 
respectively. Consequently these equations may be used to cut down the independent data 
in Taylor series coefficients of the velocity and temperature at first second and third order 
to 12, 24 and 40 components respectively. We will now redo this counting keeping track of 
the 50(3) transformation properties of all fields. 

Let us list degrees of freedom by the vector (a, b, c, d, e) where a represents the number 
of SO (3) scalars (1), b the number of 50(3) vectors (3), etc.. Working up to third order we 
encounter terms transforming in at most the 9 representation of 50(3). In this notation, 
the number of degrees of freedom in Taylor coefficients are (2, 3, 1, 0, 0), (3, 5, 3, 1, 0), and 
(4,7,5,3,1) at first, second and third order respectively. The number of equations of 
motion are (1,1,0,0,0), (2,3,1,0,0) and (3,5,3,1,0) respectively (note that the number 
of equations of motion at order n + 1 is the same as the number of variables at order n). It 
follows from subtraction that the number of unconstrained variables at zeroth, first, second 
and third order respectively can be chosen to be (1, 1, 0, 0, 0), (1, 2, 1, 0, 0), (1, 2, 2, 1, 0) and 
(1, 2, 2, 1, 1). This choice is convenient in checking the statements about the non-negativity 
of the divergence of the entropy current at third order explicitly. 



References 

[1] G. Policastro, D. T. Son, and A. O. Starinets, "The shear viscosity of strongly 
coupled N = 4 supersymmetric Yang-Mills plasma," Phys. Rev. Lett. 87 (2001) 
081601, |hep-th/0104066| . 

[2] D. T. Son and A. O. Starinets, "Minkowski-space correlators in AdS/CFT 

correspondence: Recipe and applications," JEEP 09 (2002) 042, |fiep-th/0205051 



[3] G. Policastro, D. T. Son, and A. O. Starinets, "From AdS/CFT correspondence to 
hydrodynamics," JEEP 09 (2002) 043, |aep-th/ 020505^ . 



[4] G. Policastro, D. T. Son, and A. O. Starinets, "From AdS/CFT correspondence to 
hydrodynamics. II: Sound waves," JEEP 12 (2002) 054, |hep-th/02i0220| . 

[5] C. P. Herzog, "The hydrodynamics of M-theory," JEEP 12 (2002) 026, 
|hep-th/0210126 . 



C. P. Herzog and D. T. Son, "Schwinger-Keldysh propagators from AdS/CFT 



correspondence," JEEP 03 (2003) 046, |hep-th/02i2072 



35 As T^ 1 " is not homogeneous in the derivative expansion, these equations of motion mix terms of 
different order in this expansion. 



45 



[7] C. P. Herzog, "The sound of M-theory," Phys. Rev. D68 (2003) 024013, 
[hep-th/0302086 . 



P. Kovtun, D. T. Son, and A. O. Starinets, "Holography and hydrodynamics: 
Diffusion on stretched horizons," JEEP 10 (2003) 064, |hep-th/0309213| . 

[9] A. Buchel and J. T. Liu, "Universality of the shear viscosity in supergravity," Phys. 
Rev. Lett. 93 (2004) 090602, [hep-th/03 11175 . 



[10] A. Buchel, J. T. Liu, and A. O. Starinets, "Coupling constant dependence of the 
shear viscosity in N=4 supersymmetric Yang-Mills theory," Nucl. Phys. B707 (2005) 



56-68, |hep-th/0406264 . 



[11] A. Buchel, "On universality of stress-energy tensor correlation functions in 
supergravity," Phys. Lett. B609 (2005) 392-401, [Hep-th/0408095l . 



[12] P. Kovtun, D. T. Son, and A. O. Starinets, "Viscosity in strongly interacting 

quantum field theories from black hole physics," Phys. Rev. Lett. 94 (2005) 111601, 
hep-th/0405231 . 



[13] P. K. Kovtun and A. O. Starinets, "Quasinormal modes and holography," Phys. Rev. 
D72 (2005) 086009, |hep-th/0506T8l . 



[14] P. Benincasa, A. Buchel, and A. O. Starinets, "Sound waves in strongly coupled 
non-conformal gauge theory plasma," Nucl. Phys. B733 (2006) 160-187, 



hep-th/0507026 



[15] K. Maeda, M. Natsuume, and T. Okamura, "Viscosity of gauge theory plasma with 
a chemical potential from AdS/CFT," Phys. Rev. D73 (2006) 066013, 
hep-th/0602010| . 



[16] J. Mas, "Shear viscosity from R-charged AdS black holes," JEEP 03 (2006) 016, 
[hep-th/0601144 . 



[17] O. Saremi, "The viscosity bound conjecture and hydrodynamics of M2-brane theory 
at finite chemical potential," JEEP 10 (2006) 083, |hep-th/060lT5S . 



[18] D. T. Son and A. O. Starinets, "Hydrodynamics of R-charged black holes," JEEP 
03 (2006) 052, |hep-th/060iT57| . 



[19] P. Benincasa, A. Buchel, and R. Naryshkin, "The shear viscosity of gauge theory 
plasma with chemical potentials," Phys. Lett. B645 (2007) 309-313, 
[hep-th/06 10145 . 



46 



[20] R. A. Janik and R. Peschanski, "Asymptotic perfect fluid dynamics as a consequence 
of AdS/CFT," Phys. Rev. D73 (2006) 045013, |hep-th/0512"T62| . 

[21] R. A. Janik and R. Peschanski, "Gauge / gravity duality and thermalization of a 
boost- invariant perfect fluid," Phys. Rev. D74 (2006) 046007, |hep-th/060614S. 



[22] S. Nakamura and S.-J. Sin, "A holographic dual of hydrodynamics," JHEP 09 
(2006) 020, |hep-th/0607T23| . 



[23] S. Bhattacharyya, S. Lahiri, R. Loganayagam, and S. Minwalla, "Large rotating AdS 



black holes from fluid mechanics," |arXiv : 0708 . 1770 [hep-th] 



[24] S.-J. Sin, S. Nakamura, and S. P. Kim, "Elliptic flow, Kasner universe and 



holographic dual of RHIC fireball," JHEP 12 (2006) 075, |hep-th/06 10113 . 



[25] R. A. Janik, "Viscous plasma evolution from gravity using AdS/CFT," Phys. Rev. 
Lett. 98 (2007) 022302, [hep-th/0610144 . 



[26] J. J. Friess, S. S. Gubser, G. Michalogiorgakis, and S. S. Pufu, "Expanding plasmas 
and quasinormal modes of anti-de Sitter black holes," JHEP 04 (2007) 080, 
hep-th/ 0611005| . 



[27] K. Kajantie and T. Tahkokallio, "Spherically expanding matter in AdS/CFT," Phys. 
Rev. D75 (2007) 066003, |hep-th/061222^ . 



[28] M. P. Heller and R. A. Janik, "Viscous hydrodynamics relaxation time from 
AdS/CFT," Phys. Rev. D76 (2007) 025027, [hep-th/0703243 . 



[29] K. Kajantie, J. Louko, and T. Tahkokallio, "Gravity dual of 1+1 dimensional 

Bjorken expansion," Phys. Rev. D76 (2007) 106006, |arXiv : 0705 . 1791 [hep^th] 



[30] P. M. Chesler and L. G. Yaffe, "The stress-energy tensor of a quark moving through 
a strongly-coupled N=4 supersymmetric Yang-Mills plasma: comparing 
hydrodynamics and AdS/CFT," |arXiv: 0712 .0050 [heFthJ . 



[31] P. Benincasa, A. Buchel, M. P. Heller, and R. A. Janik, "On the supergravity 

description of boost invariant conformal plasma at strong coupling," Phys. Rev. D77 
(2008) 046006, |arXiv : 0712 . 2025 Lhep^thj] . 



[32] R. Baier, P. Romatschke, D. T. Son, A. O. Starinets, and M. A. Stephanov, 
"Relativistic viscous hydrodynamics, conformal invariance, and holography," 
arXiv: 0712 .2451 |hep-th]| . 



17 



[33] S. Bhattacharyya, V. E. Hubeny, S. Minwalla, and M. Rangamani, "Nonlinear Fluid 
Dynamics from Gravity," JEEP 02 (2008) 045, |arXiv: 0712 .2456 Lhep^thl . 



[34] M. Natsuume and T. Okamura, "Causal hydrodynamics of gauge theory plasmas 
from AdS/CFT duality," |arXiv: 0712 . 2916 LhejTthl] . 



[35] M. Natsuume and T. Okamura, "A note on causal hydrodynamics for M-theory 



branes," |arXiv : 0801 . 1797 [hep-th] 



[36] A. Buchel, "Shear viscosity of boost invariant plasma at finite coupling," 
larXiv: 0801 .4421 [hep-th]| . 

[37] M. Van Raamsdonk, "Black Hole Dynamics From Atmospheric Science," 
larXiv: 0802 .3224 [hep-th] . 



[38] S. S. Gubser and A. Yarom, "Linearized hydrodynamics from probe-sources in the 
gauge- string duality," |arXiv: 0803 .0081 Lhep-thJl 

[39] L. Fidkowski, V. Hubeny, M. Kleban, and S. Shenker, "The black hole singularity in 



AdS/CFT," JEEP 02 (2004) 014, [hep-th/0306170 



[40] R. Loganayagam, "Entropy Current in Conformal Hydrodynamics," 
larXiv: 0801 .3701 |hep-th]| . 

[41] S. Hawking and G. Ellis, "The large structure of space-time," Cambridge 
Monographs on Mathematical Physics (1973). 

[42] R. Penrose, "Structure of space-time," Battelle Recontres (1967) 121-235. 

[43] V. Iyer and R. M. Wald, "Some properties of Noether charge and a proposal for 
dynamical black hole entropy," Phys. Rev. D50 (1994) 846-864, |gr-qc/9403'02£ . 



[44] I. Booth, "Black hole boundaries," Can. J. Phys. 83 (2005) 1073-1099, 
|gr-qc/0508107 . 



[45] A. Ashtekar and B. Krishnan, "Isolated and dynamical horizons and their 
applications," Living Rev. Rel. 7 (2004) 10, |gr-qc/0407042 . 



[46] E. Gourgoulhon and J. L. Jaramillo, "Area evolution, bulk viscosity and entropy 
principles for dynamical horizons," Phys. Rev. D74 (2006) 087502, |gr-qc/060705C| . 

[47] R. M. Wald, "Black hole entropy is the Noether charge," Phys. Rev. D48 (1993) 
3427-3431, |gr-qc/9307038| 



18 



[48 
[49 
[50 
[51 
[52 
[53 
[54 

[55 



[56 



[57 
[58 



V. Hubeny, A. Maloney, and M. Rangamani, "String-corrected black holes," JHEP 
05 (2005) 035, jhep-th/04ll272 . 



A. Borde, "Geodesic focusing, energy conditions and singularities," Class. Quant. 
Grav. 4 (1987) 343-356. 

T. Jacobson, G. Kang, and R. C. Myers, "Black hole entropy in higher curvature 
gravity," |gr-qc/950200S , 



B. Zwiebach, "Curvature Squared Terms and String Theories," Phys. Lett. B156 
(1985) 315. 

M. Brigante, H. Liu, R. C. Myers, S. Shenker, and S. Yaida, "The Viscosity Bound 
and Causality Violation," |arXiv : 0802 . 3318 [hep-th] . 



M. Brigante, H. Liu, R. C. Myers, S. Shenker, and S. Yaida, "Viscosity Bound 
Violation in Higher Derivative Gravity," |arXiv : 0712 . 0805 [hep-th]| . 

O. Aharony, S. Minwalla, and T. Wiseman, "Plasma-balls in large N gauge theories 
and localized black holes," Class. Quant. Grav. 23 (2006) 2171-2210, 
hep-th/0507219. 



S. Lahiri and S. Minwalla, "Plasmarings as dual black rings," |arXiv : 0705 . 3404 
Lhep-thJ 



H. K. Kunduri, J. Lucietti, and H. S. Reall, "Near-horizon symmetries of extremal 



black holes," Class. Quant. Grav. 24 (2007) 4169-4190, |arXiv : 0705~4214 



[hep-th] 



R. M. Wald, "Black holes and thermodynamics," |gr-qc/9702022 



R. Wald, "General relativity," Chicago, University of Chicago Press, 1984, 504 V- 
(1984). 

T. Jacobson, G. Kang, and R. C. Myers, "On black hole entropy," Phys. Rev. D49 
(1994) 6587-6598, |gr-qc/9312023| . 



19 



